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ABSTRACT 

Context. The nature of embedded accretion disks around forming high-mass stars is one of the missing puzzle pieces for a general 
understanding of the formation of the most massive and luminous stars. 

Aims. We want to dissect the small-scale structure of the dust continuum and kinematic gas emission toward two of the most prominent 
high-mass disk candidates. 

Methods. Using the Plateau de Bure Interferometer at ~ 1.36 mm wavelengths in its most extended configuration we probe the dust 
and gas emission at ~0.3", corresponding to linear resolution elements of ~800 AU. 

Results. Even at that high spatial resolution NGC7538IRS1 remains a single compact and massive gas core with extraordinarily high 
column densities, corresponding to visual extinctions on the order of 10 5 mag, and average densities within the central 2000 AU of 
~ 2.1 x 10 9 cnT 3 that have not been measured before. We identify a velocity gradient across in northeast-southwest direction that 
is consistent with the mid-infrared emission, but we do not find a gradient that corresponds to the proposed CH 3 OH maser disk. 
The spectral line data toward NGC7538IRS1 reveal strong blue- and red-shifted absorption toward the mm continuum peak position. 
While the blue-shifted absorption is consistent with an outflow along the line of sight, the red-shifted absorption allows us to estimate 
high infall rates on the order of 10~ 2 M Q yr~' . Although we cannot prove that the gas will be accreted in the end, the data are consistent 
with ongoing star formation activity in a scaled-up low-mass star formation scenario. Compared to that, NGC7538S fragments in a 
hierarchical fashion into several sub-sources. While the kinematics of the main mm peak are dominated by the accompanying jet, we 
find rotational signatures from a secondary peak. Furthermore, strong spectral line differences exist between the sub-sources which is 
indicative of different evolutionary stages within the same large-scale gas clump. 

Conclusions. NGC7538IRS1 is one of the most extreme high-mass disk candidates known today. The large concentration of mass into 
a small area combined with the high infall rates are unusual and likely allow continued accretion. While the absorption is interesting 
for the infall studies, higher-excited lines that do not suffer from the absorption are needed to better study the disk kinematics. In 
contrast to that, NGC7538S appears as a more typical high-mass star formation region that fragments into several sources. Many of 
them will form low- to intermediate-mass stars. The strongest mm continuum peak is likely capable to form a high-mass star, however, 
likely of lower mass than NGC75 3 8IRS1. 
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1. Introduction 

The characterization of accretion disks around young high-mass 
protostars is one of the main unsolved questions in massive star 
formation research (e.g., Beuther et al.| 2007a 2009 Cesaroni 
|et al.||200"7] |Kraus et aT1|2010) . The controversy arises around 
the difficulty to accumulate mass onto a massive protostar when 
it gets larger than 8M Q because the radiation pressure of the 
growing protostar may be strong enough to revert the gas in- 
flow in spherical accretion scenarios (e.g., Kahn 1974] |Wolfire| 
|& Ca ssinelli 1987). Different ways to circumvent this problem 
are proposed, the main two are (a) scaled-up disk accretion (e.g., 
Yorke & Sonnhalter||2002} |Krumholz et al.||2009[ |Kuiper et al.| 
2010) partially requiring initial turbulent gas and dust cores (e.g., 
McKee & Tan 2003[ ) and including ionization radiation (e.g., 
Keto|2002| l and magnetic fields (e.g., |Peters et al.|201 l| l, and (b) 
competitive accretion and potential (proto)stellar mergers at the 
dense centers of evolving massive (proto)clusters (e.g. 
|et al.|2004| [20071 |Bally & Zinnecker|2005] l. 



Bonnell 



* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). The data are available in electronic form 
at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or 
via http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A/ 



Over recent years, much indirect evidence has been accu- 
mulated that massive accretion disks do exist. The main argu- 
ment stems from massive molecular outflow observations that 
identify collimated and energetic outflows from high-mass pro- 
tostars, resembling the properties of known low-mass star for- 



mation sites (e.g., [Henning et al. 2000 Beuther et al.|2002 


Wu 


et al. 


2004 


Zhang et al. 2005 ; Arce et al.|2007 Lopez-Sepu 


Icre 


et al. 


2009 


). Such collimated jet-like outflow structures are only 



explainable with an underlying massive accretion disk driving 
the outflows via magneto-centrifugal acceleration. From a mod- 
eling approach, numeric simulations and analytic calculations 
of massive collapsing gas cores result in the formation of mas- 
sive accretion disks as well (|Yorke & Sonnhalter 2002 Kratter 

2009t |Kuiper et a 



& Matznerl 



Peters et al. 



2006 1 Krumholz et al. 



2011). Although alternative formation scenarios are 



proposed, there is a growing consensus in the massive star for- 
mation community that accretion disks should also exist in high- 
mass star formation. However, it is still poorly known whether 
such massive disks are similar to their low-mass counterparts, 
hence dominated by the central protostar and in Keplerian rota- 
tion, or whether they are perhaps self-gravitating non-Keplerian 
entities. 

While indirect evidence for massive disks is steadily increas- 
ing, direct observational studies are largely missing. This dis- 
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crepancy can mainly be attributed to the clustered mode of mas- 
sive star formation, the typically large distances and the high 
extinction. Hence, spatially disentangling such structures is a 
difficult task. While several disk candidates exist around early 
B-stars (e.g., |Cesaroni et al.||1997| [2005] |Schreyer et al l|2002 



Shep herd et al.|2001||Zhange t al. 200 21 |Chini et al 12004] |Kraus 
et al.||2010||Keto & Zhang||2010| |Fallscheer et al.|f20TT| , more 
massive O-star like systems rather show larger-scale toroid-like 
structures not consistent with classical Keplerian accretion disks 



(e.g.,|Beltran et al.|2004 
Beuther & Walsh 2008 



Beltranetal. 2011 


Beuther et al. 


2005 


Sollins et al.|2005 


Keto & Wood 


2006 



Cesaroni e t al. 2007; Fallsch eer et al.|20 09). However, the non- 
detection of Keplerian structures around O-stars does not imply 
that they do not exist, it rather indicates that they are likely on 
smaller spatial scales hidden by the toroids. Therefore, penetrat- 
ing more deeply into the central structures at the highest possible 
spatial resolution is the next step to go. 



The high-mass accretion disk candidates NGC7538IRS1 and 
NGC7538S: The source selection for such a project is driven 
by the scientific goals and the technical feasibility. The two disk 
candidates NGC7538IRS1 and NGC7538S combine the best 
of both worlds: On the one hand, they are two already well- 
studied massive accretion disk candidates in different evolution- 
ary stages at a still modest distanc e of ~2.7kpc (e.g., [Sa ndell 
et al.|2003[|Sandell & Wright|2010[|Pestalozzi et al.|2004||2006| 
Moscadell i et al.|200 9; Puga et al. 2010|), and they are easily ob- 



servable in a mosaic mode since their spatial separation is only 
of order ~80". On the other hand, with a R.A. of 23 hours and a 
Decl. of 61 degrees they are ideal Plateau de Bure Interferometer 
(PdBI) targets to be observed in long tracks resulting in the best 
achievable synthesized beams. Figure [T] presents a large-scale 
1.2 mm continuum map on top of near-infrared K-band data 
([Sandell & Sieversj [2004} |Puga et al.j [20T0] > (see also |Reid & 
|Wilson|2005D ~ 

NGC7538IRS1 has been extensively studied, and the cen- 
tral source is estima ted to a mass of ~30M o and a luminosity of 
~ 8 x 10 4 L Q (e g.,|Wmner|l976HCampbell & Thompson|l984| 
|Pestalozzi et al.|2004| >. While |Campbell| ( |1984| i and |Sandell et ak 



(20091 report an ionized jet in north-south direction, |Minier et al 
feoool and|Pestalozzi et al.|(|2004][2009]l present CH 3 OH maser 



observations indicative of an accretion disk almost perpendic- 
ular to the outflow. Partly different interpretations arise from 
mid-infrared continuum imaging (De Buizer & Minier 2005): 
They detected elongated mid-infrared emission in northwest- 
southeast direction aligned with a bipolar outflow reported in 
CO by |Davis et aE] ( |1998[ ), and even earlier in NH 3 by |Keto 



( |1991| l who also performed radiation transfer calculations for 
the region. This large-scale mid-infrared emission may stem 
from the inner walls of the outflow. The jet and outflow emis- 
sion is interpreted from speckle data as due to a precessing jet 
( |Kraus et al.| |2006j ). On smaller scales, the mid-infrared emis- 
sion appears elongated almost perpendicular to that outflow axis 
which |De Buizer & Minier| ( |2005| l interprete as an inner ac- 
cretion disk of approximate size of ~ 900AU. IKlaassen et al. 



(2009 ) also identify a velocity gradient in dense gas in northeast- 
southwest direction. Furthermore, the rare H2CO maser emis- 



sion from this source fForster et al. 1985 Pratap et al. 1992 
Hoffman et al. 2003[l is also consistent with a very young disk 



candidate. |Surcis et al.| ( |201 1) present new CH3OH and H2O 
maser observations, and their data are also consistent with a 
rotating structure in northeast-southwest direction and an out- 
flow opposite to that. Hutawarak orn & Cohen| p003 ) show the 



OH maser emission in this region. Figure [2] sketches the differ- 
ent axis and other features found in the literature. While maser, 
ionized gas and warm dust are well studied for this source, a 
good characterization of the dust and thermal gas emission was 
lacking. Recently, |Qiu et alT ( 20 1 1 1 observed the region with the 
Submillimeter Array (SMA) at 3" X 2" resolution in mm contin- 
uum and line emission, and they revealed 9 mm sources within a 
projected area of 0.35 pc. Compared to the blue-shifted absorp- 
tion observed by Keto ( 1991 1 and |Zheng et al.| ( [2001| ) in the low- 
excitation NH3 lines that is indicative of outflowing gas motions, 
Qiu et al. ( 201 l| l detected first red-shifted absorption toward the 



main mm core which they interprete as signature of ongoing in- 
fall. Furthermore, the general structure of their proposed multi- 
ple outflows is consistent with the northwest-southeast outflow 
previously reported by, e.g., Keto (1991 1; Davis et al. (1998); 



|Klaassen et al.| ( [201 1[ ). Our data now resolve the region at again 
an order of magnitude higher spatial resolution, allowing us to 
study the central core in unprecedented detail. 
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Fig. 1. Overview of the NGC7538 complex. The grey-scale 
presents the K-band image from|Puga et al.|(|2010[), and the con- 



tours show the single-dish 1.2 mm continuum data from |Sandell 
|& Sievers| ( |2004| ). The colored stars and circles show the central 
positions and FWHM of the primary beam of the PdBI at 1 .3 mm 
wavelength for NGC7538IRS1 and NGC7538S in the north and 
south, respectively. The contour levels are from 250mJybeanT 1 
to 5.25 JybeairT 1 in steps of 500 mJy beam -1 . The beam of the 
30 m observations and a scale-bar are shown at the bottom-left 
and right, respectively. 

NGC7538S is supposed to be younger than NGC7538IRS1 
but also hosts CH3OH Class n, H2O and OH maser emission 
( |Kameya et aLl[T990l [Argon et al.||2000]|Pestalozzi et al.||2006) . 
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Fig. 2. PdBI 1 .36 mm continuum images toward NGC7538IRS 1 and NGC7538S in the left and right panel, respectively. The contour 
levels start at 4<x values and continue in 8cr and 4<x steps for NGC7538IRS1 and NGC7538S (lcr values a re 29 mJy beam" 1 and 
0.7 mJy beam -1 , respectively). Several potential disk and outflow axis reported in the literature are presented ( |Davis et al.|l998 De 
Buizer & Minier|2005[|Sandell et al.|2009)|Pestalozzi et al.|2004| [20091 jSandell & Wright|2010| see Introduction for more d etails). 
The open stars, square, triangle and six-pointed star mark the positions of the OH, H2O, class 11 CH3OH and H2CO masers (Argon 
et al.!, 2000; Kamey aet al.| [T990; Pestalozzi et al.||20 06; Hoff man" et al.| 2003| >. A scale-bar and the synthesized beam are shown in 
both panels. The box zooms into the region around mml in more detail. The coordinates are relative to the phase centers given in 
sec. [2] 



About ~ 80" south of NGC7538IRS1 (Fig.[T}, it roughly coin- 
cides with a far-infrared source with luminosity ~ 1.5 X 10 4 L G 



( Werner et al.| |1979| Thronson & Harper] |1979[l c orrespond- 
ing to an early B-star. Recently, [Sandell & Wright (2010) and 
Wright et al] (|20T2|> also report the detection of NGC7538S in 



deep mid-infrared IRAC Spitzer observations at wavelengths be- 
tween 4.5 and 8 //m, as well as a weak cm continuum source 
likely stemming from a thermal jet. Sandell et al. (2003 ) resolved 
a 30000 AU rotating structure with a bipolar outflow emanating 
perpendicular to that by means of BIMA interferometric mm ob- 
servations . Wright et al. ( 2012 1 resolved that structure into three 



separate mm sources. Recent interferometric observations (be- 
tween ~ 3" and 8" resolution) in different molecular line tracers 
largely confirm this picture (Sandell & Wright 2010). However, 
their work indicated that many of the molecules are affected by 
the jet/outflow, and clear rotational signatures were hard to iso- 
late. While the previous data are consistent with a rotating struc- 
ture, the spatial resolution was not sufficient to analyze the disk 
candidate in detail. Figure [2] again sketches the main features 
found in the literature. 

The different ages of the two targets make them ideal candi- 
dates to investigate also disk evolutionary properties within the 
same observations. 



2. Observations 

The two sources were observed in two tracks - A and B con- 
figuration on January 26th, 2011, and February 10th, 2011, 
respectively - in two fields where each field was centered 
on NGC7538IRS1 and NGC7538S (see Fig. [TJ. The phase 
centers for the two fields were R.A (J2000.) 23h 13m45.36s, 
Dec. (J2000.0) 61° 28' 10.55" and R.A. (J2000.) 23h 13m44.86s 
and Dec. (J2000.0) 61° 26' 48.10". Phase calibration was con- 



ducted with regularly interleaved observations of the quasars 
2146+608, 0059+581 and 0016+731. The bandpass and flux 
were calibrated with observations of 3C345, 3C273 and 
MWC349. The absolute flux level is estimated to 20% accuracy. 

The continuum emission was extracted from broad band 
data obtained with the WIDEX correlator with four units and 
two polarizations covering the frequency range from 217.167 
to 220.836 GHz. NGC7538IRS1 is such an extremely line- 
rich source that barely any line-free region exists in the spec- 
trum. Therefore, extracting a line-free continuum is difficult, 
and we produced our continuum map from the whole bandpass. 
Although there obviously is some line contamination in the final 
continuum image, the continuum level itself (> 1.8 Jy beam -1 ) 
is so extraordinarily high that the relative contribution from the 
lines is negligible. To get a quantitative estimate of the line con- 
tamination, we also produced a continuum image from only a 
very small line-free region (~ 71MHz) between the CH3CN 
lines. The peak flux in this image is only 1.1% below the peak 
flux in the continuum map produced from the whole band- 
pass. Obviously the line contamination is negligible and we use 
the map based on the full bandpass because of the lower rms 
noise. NGC7538S is less line-rich, and therefore, we could pro- 
duce the continuum image from the line-free parts of the spec- 
trum. A comparison of the continuum images for NGC7538S 
with and without line contamination shows that they are almost 
identical. Although NGC7538IRS1 exhibits more line emis- 
sion, the previous comparison also indicates that including the 
lines in the NGC7538IRS 1 continuum bandpath only marginally 
changes the real fluxes. The full bandpath spectral line data with 
a chemical analysis will be presented in a forthcoming paper. 
The lcr continuum rms for NGC7538IRS1 and NGC7538S are 
29mJybearrT 1 and 0.7 mJy beam" 1 . The difference in rms can 
be explained by the fact that for sparse antenna interferometers 
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like the PdBI, the rms is usually not the thermal rms but it is 
dominated by the side-lobes of the strongest source in the field. 
And since NGC7538IRS1 is far brighter than NGC7538S, also 
the rms for NGC7538I is significantly higher. 



Table 1. Observed spectral lines 



CH 3 CN(12 k -1 1 k ) toward continuum peak (offset 0.1 "/-0.1 5") 
I 1 1 1 



Freq. 




Mol. 


EJk 


(GHz) 






(K) 


218.222 


H 2 CO(3„,3 


-2 ,i) 


21 


218.298 


HCOOCH 3 (17 3 14 - 


I6343) 


100 


218.325 


HC 3 N(24 - 23) 


131 


218.440 


CH 3 OH(4 22 


-3 U2 ) 


46 


218.460 


NH 2 CHO(10, 9 


-9i, 8 ) 


61 


218.476 


H 2 CO(3 2 , 


-2 2 ,,) 


68 


218.903 


OCS(18- 17) 


100 


220.167 


HCOOCH 3 (17 4>13 - 


164,12) 


103 


220.178 


CH 2 CO(llni - 


IO1.10) 


77 


220.190 


HCOOCH 3 (17 4 ,i 3 - 


164,12) 


103 


220.594 


CH 3 CN(12 6 


-He) 


326 


220.641 


CH 3 CN(12 5 


-Us) 


248 


220.679 


CH 3 CN(12 4 


-H4) 


183 


220.709 


CH 3 CN(12 3 


-lis) 


133 


220.730 


CH 3 CN(12 2 


-H2) 


98 


220.743 


CH 3 CN(12! 


- Hi) 


76 


220.747 


CH 3 CN(12 


-Ho) 


69 



To extract kinematic information, we put several high- 
spectral resolution units with a nominal resolution of 0.3 12 MHz 
or 0.42 km s _I into the bandpass covering the spectral lines listed 
in Table [l] The spectral line rms for 0.5kms~' wide spectral 
channels measured in emission-free channels is SmJybeairr 1 
and 7mJybeairT 1 for NGC7538IRS1 and NGC7538S, re- 
spectively. The v lsr for NGC7538IRS1 and NGC7538S are ~ 

-57.3 km s" 1 and 56.4kms— 1, respectively (Gerneret al. in 

prep., |van der Tak et al.|200"0l|Sandell & Wright|2010| i. The data 
were inverted with a "robust" weighting scheme and cleaned 
with the dark algorithm. The synthesized beam of the final con- 
tinuum and line data is ~ 0.31" x 0.29" (PA 1 10°). 



3. Results 

3.1. Continuum emission 

Zooming in from the large-scale emission (Fig. [TJ, Figure [2] 
presents the small-scale structure of the region in the 1.36 mm 
continuum emission at a spatial resolution of ~ 0.3" or ~ 
800 AU. While NGC7538IRS1 remains a single source with 
peak flux in excess of 1.8 Jy beam -1 , NGC7538S is resolved 
into several sub-sources labeled mml to mm3. NGC7538S mml 
shows even additional fragmented substructure which we la- 
bel as mmla and mmlb. The sub-source mmla is elongated 
in approximately the north-south direction, and it is likely that 
this elongation corresponds to an unresolved substructure again. 
Comparing the high-resolution mm continuum data with the 
near-infrared image by Puga et al. (2010 1, the mm peak in 
NGC7538IRS1 is clearly associated with the main infrared 
source IRS1, whereas the three mm peaks in NGC7538S have 
no near-infrared counterpart. However, NGC7538S mml has re- 
cently been detected by Spitzer at wavelengths between 4.5 and 
8//m USandell & Wright) |20T0l [Wright et al.[ |2012) It should 



also be noted that the 8 additional sources reported by Qiu et al. 
( 201 l| l within an area of 0.35 pc are not detected by our higher 
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Fig. 3. CH 3 CN(12* - ll k ) spectra for k = 0...3 toward the 
1.36 mm continuum peak (top panel) and a reference position 
offset by 0.2"/ - 1.1" from NGC7538IRS1. The velocity of rest 
is always marked with a vertical line. 



to our smaller primary beam (FWHM of ~ 22") as well as the 
lower brightness sensitivity one automatically achieves when 
going to higher spatial resolution (our 3cr continuum rms of 
87mJybearrT 1 corresponds to an approximate brightness sen- 
sitivity of -24 K). 

Table [2] presents the measured peak and integrated fluxes 
of the sub-sources shown in Fig. [2] The integrated fluxes 
are measured within the 4<x contours. The single-dish 1.2 mm 
MAMBO data shown in Figure [T] exhibit peak fluxes of 5702 
and 2872 mJy beam 1 for NGC7538IRS1 and NGC7538S, re- 
spectively. Comparing these numbers to the integrated fluxes 
we measure with the PdBI (Table [2}, we find that toward 
NGC7538IRS1 only ~ 48% of the flux is filtered out with the 
interferometer. |Qiu et al7| ( |201 l| l also measure with the SMA at 
about an order of magnitude lower spatial resolution (3" x 2") 
an integrated flux of 3.6 Jy, only ~20% higher than our fluxes 
measured with ~ 0.3" resolution. In comparison to that, to- 
ward NGC7538S approximately 90% of the single-dish flux is 
missing in the PdBI data. This large difference indicates that 
NGC7538IRS1 is extremely concentrated toward the central mm 
continuum peak whereas NGC7538S exhibits emission on much 
larger scales. Since NGC7538IRS1 has significant amounts of 
free-free emission (e.g., Pratap et al. 1992 Keto et al. [2008 



Sand ell et al. [20091 ), we correct the fluxes for that contribution 
in Table |2j As shown by |Keto et ak] ( |2008| l, several Hn region 
models can fit the data, and the exact free-free contribution is 
hard to isolate. Here, we assume ~1000mJy to be produced by 
the free-free emission, the rest is attributed to the dust emis- 



resolution PdBI observations. This differences can be attributed 



sion. Recently, Wright et al. ( 2012 1 report the cm free-free fluxes 
from NGC7538S mmla, and following their approach we cor- 
rect 8 mJy free-free flux contribution in Table |2]as well. 

Assuming optically thin emission from dust following the 
standard approach by Hildebrand (1983) we can estimate 
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Fig. 4. The left and right panel present in color the 1st and 2nd moment maps in CH3CN(126 - lU) toward NGC7538IRS1. The 
contours show the 1 .36 mm continuum emission with the same contour levels as in Fig. [2] A scale-bar and the synthesized beam are 
presented in each panel. The line in the left panel outlines the axis for the pv-cuts in Figures [5] and [6] 



Table 2. Millimeter continuum properties 



Source 


S peak 
mJy 
beam -1 


int 

mJy 


M 
M G 
H83 


N 
io 25 

cm -2 
a 


M N 

w cm 1 

OH94" 


7538IRSlmml 


1861 


2990 










7538IRSlmml' v 


861 


1990 


115 


18 


43 


7 


7538Smml 


24.1 


132 










7538Smml w 


16.1 


124 


38 


1.8 


14 


0.7 


7538Smmla 


24.1 


95 










7538Smmla M 


16.1 


87 


27 


1.8 


10 


0.7 


7538Smmlb rf 


18.0 


37 


11 


2.0 


4 


0.7 


7538Smm2' ; 


20.5 


76 


23 


2.3 


9 


0.9 


7538Smm3 rf 


21.4 


67 


21 


2.4 


8 


0.9 



very high, similar values have been reported in the past at cor- 
respondingly high spatial resolution (e.g., Beut her et al.|2 007c 



n H83: Hildebrand ( 1983 ); OH94: |Ossenkopf & Henni ngl 1994) 

* Corrected for free-free flux contribution. 

c Assumed temperature 245 K. 

d Assumed temperature 50 K. 

e Measured within the 4cr contours. 



gas masses and column densities at an assumed temperature. 
Because NGC7538IRS1 is a strong infrared source and hot core, 
following |Qiu et aL] ( |201 1[ ) we assume a dust temperature of 
245 K for that source. NGC7538S is supposedly younger and 
colder, and we assume a dust temperature of 50 K for the cor- 
responding sub-sources. Regarding the dust properties, we cal- 
culate the masses and column densities following Hildebrand 
( fT98"3] l (H83) on the one hand, and |Ossenkopf & Henning| ( |1994| l 
(OH94) for thin ice mantles at densities of 10° cm" 3 on the other 



hand. The gas-to-dust ratio is taken as 186 following Draine et al. 
p007l > and |Jenkins| ( p004l i. 

Depending on the dust properties, toward NGC7538S we 
find core masses between 4 and 38 M Q and column densities be- 
tween 0.7 x 10 25 and 1.8 X 10 25 cm , corresponding to visual 
extinctions on the order of 10 4 mag. While such extinctions are 



Rodon et al. 2008 ). One should keep in mind that such high ex- 
tinction values are only found at the highest spatial resolution 
achievable with interferometers. At lower resolution, the emis- 
sion smears out and lower values are found. Regarding the core 
masses in NGC7538S, at first sight they do not appear extraordi- 
narily high, however, considering that approximately 90% of the 
gas are filtered out on larger scales, we only observe the densest 
structure that is embedded in a much larger gas reservoir. 



The situation is considerably different for NGC7538IRS1 
where excessively high column densities on the order of 
10 26 crrT 2 are found (corresponding to visual extinctions above 
10 5 mag), as well as core masses between 43 and 115M 
(depending on the dust properties) within a projected size of 
~2000AU. To the authors' knowledge this is an extraordinary 
concentration of mass within small spatial scales and will be dis- 
cussed in more detail in section 14.1.1 1 



For comparison, we can also calculate the total gas masses 
of the two regions based on the single-dish data. As approx- 
imate clump sizes, we integrate the flux in the area within 
the 750 mJy beam contour in Figure [T] For NGC7538IRS1 and 
NGC7538S, we get integrated 1.2 mm fluxes of 20.9 and 6.2 Jy, 
respectively. Following the same approach as above, we can cal- 
culate the gas masses for the two dust models H83 and OH94. On 



these large scale we use the temperature estimates from Sandell 
|& Sievers| ( p004l > who estimate 75 and 35 for NGC7538IRS1 
and NGC7538S. With these numbers we get total gas masses for 
NGC7538IRS1 and NGC7538S of 2512 and 1757 M (H83) or 
1011 and 706 M G (OH94), respectively. 
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3.2. Spectral line emission 

All spectral lines listed in Table [T were detected toward 
NGC7538IRS1 and most of them also toward NGC7538S. 
While also velocity gradients are identified in both regions, the 
detailed spectral line signatures between NGC7538IRS1 and 
NGC7538S are considerably different, in particular we detect 
strong absorption signatures toward NGC7538IRS1 but not to- 
ward NGC7538S. 



3.2.1. NGC7538IRS1 

Figure [3] presents the CH 3 CN(12 A - 11^) (0 < k < 3) spec- 
tra toward the mm continuum peak as well as toward a position 
approximately 1.1" offset south of the continuum peak. While 
the offset spectrum is a typical CH3CN emission spectrum, the 
spectrum toward the continuum peak is dominated by absorp- 
tion features. While absorption features in interferometric data 
should always be taken with caution because missing flux prob- 
lems can also artificially produce such features, the fact that we 
see the absorption only toward the continuum peak but not to- 



ward an offset of only ~ 1.1" is a strong indicator for the ab- 
sorption being a real feature. While [Qiu et al. (201 II) reported 



only redshifted absorption in dense gas tracers at a spatial res- 
olution of ~ 3" x 2", and [Ketol ( [T9^] ) and |Zheng eFaTlpOOT) 
found only blue-shifted absorption in lower-density NH3 lines, 
Figure [3] already shows that at the high-spatial resolution of our 
observations, the dense gas shows blue- and red-shifted compo- 
nents simultaneously. 

To identify potential rotational symmetries in these data, 
Figure [4] presents the 1st and 2nd moment maps (intensity- 
weighted peak velocities and line-widths) of our highest ex- 
cited CH 3 CN(12 6 - 11 6 ) line (EJk = 326 K). Although these 
moment maps are affected by the absorption close to the con- 
tinuum peak, both maps clearly identify a velocity gradient 
in northeast-southwest direction, consistent with the previous 
lower-resolution Submillimeter Array observations in OCS and 
( |2009) . 



SO2 by Klaassen et al 



While the 1st moment map ex- 
hibits a blue-red velocity gradient extending about lOkrns -1 
which is perpendicular to the northwest-southeast outflow struc- 
ture reported by Davis et al. ( 1998 1 and Qiu et al. ( |2011 1, it is 
interesting that also the line-width map shows a significant line- 
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width increase close to approximately this axis while the line- 
widths are considerably smaller northwest and southeast of that. 
The 2nd moment map in Figure |4] gives visually the impression 
of a disk-like structure, however, again this needs to be taken 
with caution because that signature can be affected by the ab- 
sorption of the gas against the strong continuum. 

Figures [5] and [6] now show the position-velocity diagrams 
along the northeast-southwest cut outlined in Figure [4] The typ- 
ical hot molecular core and high-density gas tracers CH3CN and 
HCOOCH3 exhibit absorption signatures that are dominated by 
a red-shifted component but show some blue-shifted absorp- 
tion as well. This signature appears rather independent of the 
excitation temperature because the shown k — 2 and k = 6 
CH3CN(12,t - 1 lk) components cover a range in excitation tem- 
peratures E u /k of 250 K (see Table [TJ. As expected, the more 
optically thin isotopologues CH^CN does not exhibit such ab- 
sorption features. While the red-shifted emission part of these 
spectra is consistent with a typical Keplerian rotation structure, 
the blue-shifted part of the emission spectrum does not show 
such a signature. Figure [5] also presents a Keplerian curve for 
a 30 M Q central object (see section [T}, again showing the rea- 
sonable agreement on the red part of the spectrum but not on 
the blue side. It is also interesting to note that the higher excited 
CH3CN(12g- 1 16) shows on average abetter agreement with the 
Keplerian curve than the lower excited CH3CN(122 — 1 12) line. 
This is likely due to the fact that the higher excited line traces gas 
closer to the star which hence exhibits higher velocities. In con- 
trast to that, one can also argue that the non-correspondence of 
the blue part of the spectrum with Keplerian rotation is not much 
of a surprise because Keplerian rotation implies that the whole 
structure is dominated by the central object. This is clearly not 
the case considering that the central star should have a mass of 
~30 M (see Introduction) and the gas mass of the central struc- 



ture traced by the mm continuum emission is of that order or 
even higher as well (Table |2j. 

The position-velocity diagrams presented in Figure [6] show 
slightly different absorption signatures. While almost all of them 
show red-shifted absorption as well (except CH2CO), the blue- 
shifted absorption is at least as strong as that, if not stronger. 
From a blue-shifted perspective, the lowest-excited H2CO line 
shows a particularly interesting feature because in addition to the 
blue-shifted component at around -59 km s -1 , it exhibits another 
blue-shifted component at even more negative velocities around 
-65 kms~'. A similar feature was recently also reported in the 1 
and 2 cm H 2 CO lines (AAS poster by |Yuan et al.|2011| >. 

A different way to investigate the various absorption features 
is presented in Figure [7] where we show the spectra of the vari- 
ous lines extracted directly toward the mm continuum peak posi- 
tion. The absorption features discussed in the paragraphs above 
are exactly recovered there. Many of the spectra clearly show a 
double-dibbed signature red- and blue-shifted around the v'i sr . To 
check whether the additional higher-velocity absorption compo- 
nent in the H2CO line is real absorption against the continuum 
or rather due to missing flux on larger scales, similarly as shown 
for CH3CN in Figure[3] we also extracted the H2CO spectrum to- 
ward the position 0.2"/- 1.1" to the south. And like for CH3CN, 
the H2CO spectrum exhibits a pure and "normal" emission spec- 
trum at that position. This implies that the additional absorption 

at 65 km s _1 should be real. Implications of the observed red- 

and blue-shifted absorp tion toward NGC7538IRS1 will be dis- 
cussed in section l4.2.1l 

3.2.2. NGC7538S 

Toward the second region NGC7538S we clearly detect all 
CH3CN lines, as well as the spectral lines from OCS, HC3N, 
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Fig. 7. Spectra toward the 1.36 mm continuum peak of NGC7538IRS1 for selected lines as marked in each panel. The velocity of 
rest is marked by a line. 



H2CO and CH3OH. In contrast to that, HCOOCH3, NH2CHO clearly detected molecules and spectral lines, it is interesting that 
and CH2CO are barely detected. There is only a tentative de- all of them are detected toward the two mm sub-peaks mml 
tection of the latter two molecules toward mm2. Regarding the and mm2 but none of them toward the third mm peak mm3. 
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This already indicates peculiar chemical and evolutionary differ- in the spectral line emission but no molecular line is found to- 
ences between mml and mm2 on the one side and mm3 on the ward mmlb. A detailed spectral and chemical analysis of all the 
other side. Furthermore, within mml we always detect mmla other broadband line data we observed simultaneously will be 
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presented in a forthcoming paper. Here we concentrate on the 
kinematics of the mm peaks mml and mm2. 

Figure [8] presents the 1st and 2nd moment maps (intensity 
weighted peak velocities and line widths) of respective lines to- 
ward that region. While the lower excited lines like F^CO^^ - 
224) or CH30H(42,2 - 3i,2) show also a bit more extended emis- 
sion, molecular emission from high-density tracers like CH3CN 
or HC3N are largely confined to mml, mm2 and their close envi- 
ronment. While there is no obvious velocity difference between 
mml and mm2, toward both sub-peaks we detect in all lines ve- 
locity gradients across the mm continuum peaks, for mml al- 
most north-south and for mm2 in northeast-southwest direction. 
Toward mml it is interesting to note that the molecular emis- 
sion to the north is confined almost to the same region as the 
mm continuum emission whereas the molecular line data ex- 
tend significantly outside the southern Act contours of the mm 
continuum emission. Since the main mmla peak is close to the 
northern edge of that structure, the velocity structure of the gas 
is rather asymmetric with respect to that peak. The most blue- 



shifted gas almost peaks toward mmla and the further one goes 
south, the more redshifted the gas gets. For the higher density 
lines of CH 3 CN(12* - 1 1*), OCS(18-17) and HC 3 N(24-23), the 
line-width or 2nd moment peak is also not exactly toward the 
main mm continuum peak but a little bit to the south, almost at 
the tip of the central, north-south elongated mm continuum con- 
tour. As discussed in section 3.1 we cannot properly resolve a 
secondary component there, nevertheless, it appears likely that 
this slightly elongated structure will resolve into a binary system 
at even higher spatial resolution. 

Figure [9] presents position velocity cuts of selected lines 
along the axis marked for mml in Figure [8] While the emission 
appears relatively symmetric around the vi sr , as already men- 
tioned above, the velocity is not symmetric around the main mm 
peak mmla which is put at offset 0. Even if one shifts the center 
by ~ 0.25" south toward the peak of the 2nd moment maps, it 
still does not appear as a symmetric position velocity cut. The 
data clearly confirm that the most blue-shifted gas is centered 
on mmla and the red-shifted emission continuously moves to 
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Fig. 10. Position velocity cuts through NGC7538S mm2 in several lines as marked in each panel along the axis shown in Fig. [8 
respectively. The Visr is marked. The full and dotted lines in the OCS panel correspond to Keplerian rotation curves around a 1 anc 
10M o star, respectively. The units of the wedges are Jy beam -1 . 



the south. Furthermore, the pv-diagrams do not exhibit any sig- 
nature of Keplerian rotation. These signatures indicate that the 
observed velocity structure from mml unlikely stems from ro- 
tation. Since the jet axis is aligned approximately in northwest- 
southeast direction (Figure |2j, not much offset from the main 
velocity gradient observed here, it may well be that the velocity 
gradient is strongly influenced by the central jet and outflow. 

In comparison to mml, Figure 10 shows the position- 
velocity cuts through mm2 along the axis shown in Figure [8] 
Since mm2 is much smaller in spatial extend and only barely 
resolved by our observations, the pv-diagrams also exhibit less 
prominent signatures of velocity gradients. Nevertheless, a ve- 
locity gradient is identifiable, most prominently in the OCS(18- 
17) line. While the structure is too small for a more detailed 
analysis, it is interesting to note that the OCS(18-17) data are at 
least consistent with Keplerian rotation around a ~ 1 M central 
object, whereas the expected rotation curve of a more massive 
object represents the data significantly worse. While one should 
not take these masses at face-value, they nevertheless indicate 



that the central mass in NGC7538S mm2 is significantly lower 
than that in NGC7538IRS 1 . 



A different way to investigate the spectral structure of the 
two sources is again via directly looking at the spectra toward 
the peak positions (Figures 11 and 12 1. While the spectra to- 
ward mm2 exhibit more or less Gaussian shapes around the vi sr , 
this is not the case for the spectra extracted toward mm la. The 
mml a spectra are strongly dominated by the blue-shifted gas 
which was already identified in pv-diagrams, but we see an ad- 
ditional red-shifted component that is separated by a little flux- 
depression (not absorption) around the vi sr . This will be dis- 
cussed in more detail in section l4.2.2l 
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Fig. 11. O^CN^ - ll fe ) spectra (0 < k < 3) toward 
NGC7538S mm la (black) and mm2 (red). 

4. Discussion 

4.1. Continuum structure 
4.1.1. NGC7538IRS1 

While fragmented sub-sources like those in NGC7538S have 
been observed regularly with interferometers toward high- 
mass star-forming regions (e.g., [Beuther et al. 2007b| [Zhang 
|et al.| 2009| ), the very massive, centrally condensed and unre- 
solved source NGC7538IRS1 is peculiar and less typical (e.g., 
Bontemps et al.|2010 1. In particular, the fact that one observes a 



near-infrared source toward that position but at the same time has 
column densities in excess of 10 26 crrT 2 (corresponding to visual 
extinctions in excess of 10 5 mag) is surprising. This seems to in- 
dicate that the outflow/jet from the source must be very close to 
the line of sight allowing the infrared radiation to escape through 
the outflow cavity. This is reminiscent to similar sources like 
W3IRS5 or G9. 62+0. 19 where also the infrared sources are de- 
tected in spite of high column densities because of the alignment 
of the outflow close to the line of sight (e.g., Rodon et al.[2 008 ; 
|Hofner et al.|200T]|Linz et al.|2005| l. 

Furthermore, it is astonishing how much mass is concen- 
trated in a very small area in NGC7538IRS1. Tableland Figure 
|2]show that at least more than 40 and potentially even more than 
100 M are concentrated within a source with projected diame- 
ter of approximately 2000 AU. To the authors knowledge, there 
is no other region known with so much mass in such a small area. 
Assuming 50 M within a sphere of diameter of ~2000 AU, this 
corresponds to a mean H2 density of ~ 2.1 x 10 9 cirT 3 which is 
also high compared to other star-forming regions. It is also about 
2 orders of magnitude larger than the average densities derived 



by |Qiu et al. ( 201 1) . This difference can largely be attributed to 
the much smaller size of the core in the PdBI data (a factor 3 
in radius) as well as different assumptions in the mass calcula- 
tions (we adopted a slightly different free-free contribution and 
gas-to-dust ratio). 

|Puga et aL (2010 1 report an age spread between 0.5 and 



2.2 Myrs for the surrounding infrared cluster whereas the vari- 
ous signatures of ongoing star formation toward the central mm 
and infrared source (e.g., maser emission, outflows) indicate that 
the central and most massive object is still in an active star for- 
mation process. It appears that in that region the most massive 
star forms last compared to the lower-mass population, similar 

( |2006| l or |Wang et al 



e.g., 



Kumar et al. 



to other studies like 

( |201 l| l. Furthermore, the strong concentration of mass within a 
single object without much further fragmentation (except those 
on larger scales as reported by Qiu et al. 201 1 see also Bontemps 
|et al.|2010] i is consistent with a scaled-up low-mass star forma- 
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Fig. 12. CH 3 CN(12 5 - 1 1 5 ) and HC 3 N(24 - 23) spectra toward 
NGC7538S mmla and mm2. 



tion scenario for the formation of high-mass stars (e.g., McKee 



& Ta n 2002 2003 |Krumholz et aL]|2007| 12009] |Kuiper et al. 
201 Op . 



4.1.2. NGC7538S 

The substructure one finds in NGC7538S is indicative of hi- 
erarchical fragmentation on different scales. While the single- 
dish mm continuum data show only one large-scale gas clump 
with a projected diameter of ~0.5pc (e.g., Figure [T] Sandell & 
|Sieve rs 2004 ; Reid & W ilson 2005), first interferometric obser- 
vations already revealed an elongated gas clump with an extend 
of ~30000 AU~0.15pc that showed velocity signatures indicat- 
ing rotation. Recently, Wri ght et al.| ( |2012] l resolved that elon- 
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gated structure into three mm continuum sources which we con- 
firm here at even higher spatial resolution. Our new data now 
also indicate that mml splits up likely in >3 sources forming a 
trapezium-like system (e.g.,[Ambartsumian 1955; Megeath et al.| 



2005 Rodon et al. 20081. The hierarchical fragmentation ob- 



served in NGC7538S resembles the structures recently discussed 
by |Zhang etlE] ( |2009] l and |Wang efaL] ( |20"TT] l for the infrared 
dark cloud G28.34, although on projected smaller spatial scales 

inNGC7538S. ' 

Following the simple toy-model outlined in Beuther et al. 
(|2012|l, in a cluster-forming scenario with a typical |Kroupa 



(2001 



initial mass function and a star formation efficiency of 
30%, one needs approximately a 1000 M Q initial gas clump 
to form a cluster with at least one 20 M Q high-mass stars. As 
estimated in section [XT] NGC7538S fulfills that criterium (and 
NGC7538IRS1 even more), and while mm2 supposedly does not 
form a high-mass star but rather a low- to intermediate-mass ob- 
ject (see section 3.2.2 1, the higher gas mass (Table |2]i and the 
other high-mass star formation indicators discussed in the intro- 
duction and outlined in Figure|2]indicate that mml likely forms 
a high-mass star within this still very young cluster-forming re- 
gion. 



4.2. Kinematics 
4.2.1. NGC7538IRS1 

The exceptional red- and blue-shifted absorption features pre- 



sented in section 3.2.1 give various insights into the physical 
processes around that source. As discussed in section 4.1.1 the 
jet/outflow from this source has to be oriented approximately 
along the line of sight. Therefore, blue-shifted absorption fea- 
tures against the continuum peak position should be associated 
with expanding motions from the jet/outflow. Also the fact that 
we observe in the line with the lowest excitation temperature 
(H2CO(3o,3 - 2o,i) with E„/k ~ 21 K, see also TablefTh an addi- 
tional absorption feature at even more blue-shifted velocities is 
consistent with this picture. The lower-excited line traces colder 
gas further away from the source, and many outflows are known 
to exhibit Hubble-like velocity structure where the velocity in- 
creases with distance from the source (e.g., Arce et al. 2007) . 



Hence colder gas further outside should show absorption at ve 
locities further blue-shifted than closer to the source. 

In contrast to that, the red-shifted absorption is indicative 
of infalling gas. Following the approach outlined in Qiu et al. 
( |201 l| l, assuming a spherical infall geometry one can estimate 
mass infall rateaH according to M m — 47rr 2 pvj n where M m and 
Vin are the infall rate and infall velocity, and r and p the core 
radius and density. The latter two values are ~ 1000 AU and 



2.1 x 10 9 cm 3 , respectively (see section 



4.1.1 



As infall ve- 
locity we use 2.3 km s which corresponds to the difference be- 
tween the peak of the red-shifted absorption at 55 km s _1 and 

-57.3 km s _1 (Gerner et al. in prep 



van der 



the vi sr at around 

|Tak et alj2 000). With these numbers, we derive an infall rate es- 
timate of M[ n ~ 7 x 10~ 2 M yr~' . This is approximately a factor 



20 larger than the rate estimated by Qiu et al. ( [2011) . This dif- 
ference is due to the higher spatial resolution of our data where 
we resolve the core on smaller scales (radius of 1000 AU here 
compared to 3000 AU in Qiu et al.|201 l| l which additionally re- 
sults in higher average densities or the central core (see section 
4.1.1 1. Considering that the accretion does not occur in a spheri- 



cal mode over 4-n but rather along a flattened disk structure with 
a solid angle of O, the actual disk infall rates Mdisk.in should scale 
like Mdisk.in = H x Mj n . Based on the simulations by Kuiper et al. 
( |2012| l and R. Kuiper (priv. comm.), such outflow covers approx- 
imately 120 degree opening angle and the disk 60 degree (to be 
doubled for the north-south symmetry). Since the opening angle 
does not scale linearly with the surface element, full integration 
results in ~50% or ~ 2n of the sphere being covered by the disk. 
This results in disk infall rates of Mdisk.in ~ 3.5 x 10~ 2 M yr~', 
still very high and in the regime of accretion rates required to 



form high-mass stars (e.g.,|Wolfire & Cassinelli 1987 McKee & 



|Tan| [2003). Although we cannot prove that the gas falls in that 
far that it can be accreted onto the star (and does not get reverted 
by the innermost radiation and outflow pressure), such high in- 
fall rates should be a pre-requisite to allow accretion even when 
the central high-mass star has ignited already (e.g., Keto 2003; 
|Kuiper et al.|20T0l|20TT i. 



An additional caveat arises from the potential contribution 
of the accretion luminosity to the total luminosity of the re- 
gion. If one used the infall rates as actual accretion rates M acc 
and estimated the accretion luminosity L acc via the classical 
iacc = CM ^ M " CL (with G the gravitational constant, M* and R t 
the estimated stellar mass of 30 M and a stellar radius follow- 
ing Hosokawa & Omukai 2009 ), one would derive unreasonably 



high accretion luminosities in excess of the measured luminos- 
ity. Therefore, some parameters in this equation have to be dif- 
ferent. Most likely this is the accretion rate because not all gas 
will fall on the star but a large fraction will likely be expelled 
again by the energetic outflow. Estimating that ratio is out of the 
scope of this paper. Nevertheless, the data indicate that a signifi- 
cant fraction of the measured luminosity may still stem from the 
accretion processes. 

As outlined in the introduction, the disk orientation in this re- 
gion has been subject to intense discussion. Although the 1st mo- 
ment map of NGC7538IRS1 is distorted from the absorption to- 
ward the peak, our data clearly support the orientation of the disk 
along a northeast-southwest orientation that was also proposed 
by|De Buizer & Minierl ( |2005| l, |Klaassen et alT] ( |2009| l or |Surcis| 
et al. ( |2011| l. We do not find signatures of rotational motion 



along the more east-west oriented structure that was proposed as 
a disk from CH3OH maser observations (Pestal ozzi et al.| |2004 
2009 ). However, we cannot exclude that the dense gas northeast- 
southwest oriented emission and the east-west CH3OH masers 



1 Since one does not know whether the gas gets actually accreted or 
not, we prefer the term "infall rate" rather than "accretion rate". 



belong to the same torus-disk structure as suggested by Surcis 
|e^(|20TTj. 



4.2.2. NGC7538S 

The spectral line signatures in NGC7538S vary considerably 
among the three main sub-sources. While mml and mm2 are 
strong line emitters, mm3 shows no line emission in any of 
the discussed lines of this project. Contrary to that discrepancy, 
the continuum emission from mm2 and mm3 is very similar in 
size, column density and mass. Therefore, it is most likely that 
the spectral line differences are real chemical differences among 
the two sub-sources mm2 and mm3. It is tempting to interprete 
that these differences are due to different evolutionary stages. 
While mm2 shows also rotational signatures and is already a 
star-forming core, mm3 may well still be in a younger and star- 
less phase. In this picture, sources that are separated by less than 
10000 AU and that are embedded within the same large-scale gas 
clump may not evolve coeval at all. What are the physical rea- 
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sons for this evolutionary differences? Unfortunately, our data 
do not allow us to draw conclusions on that point. 

The additional blue-redshifted gas components visible in the 
spectra toward mmla (Figures [TT| and [12} indicate that while 
the pv-diagrams (Figure |9]l are dominated by the larger struc- 
ture encompassing also the dust elongation toward the south, 
there exist additional velocity structure toward the peak mmla. 
Although we do not spatially resolve that substructure, it may 
well stem from a smaller embedded disk centered on mmla 
which is also the likely driving source of the jet in that region 
(Fig. [2}. Since the continuum source mml is already resolved in 
at least two sub-sources (mmla and mmlb), and mmla is elon- 
gated indicating the existence of an additional source, it is likely 
that mml hosts a multiple system where large-scale kinematic 
structures are present (best visible in the pv-diagram in Figure 
[9]) as well as potential small-scale rotational structure around in- 
dividual sub-sources only identified in the spectra (Figures 1 1 



and 12 1. Comparing the larger almost north-south velocity gra- 
dient in mml with the orientation of the jet in approximately 
northwest-southeast direction, it appears that the dense gas kine- 
matics in mml are strongly influenced by the jet. This makes the 
identification of rotational signatures even harder. It is likely that 
higher spatial resolution as well as spectral lines sensitive only to 
the innermost region around the central protostar are needed to 
disentangle the rotational structure around mmla from the kine- 
matic signatures caused by the jet. 



Going to larger spatial scales, Sandell et al. ( 2003| l al- 
ready identified a velocity gradient across the whole 30000 AU 
structure that encompass the three regions mml to mm3 ap- 
proximately along the connecting axis of the sub-sources. 
Interestingly, Sandell et al. ( 2003[ ) find that the large-scale ro- 
tation is consistent with Keplerian motion. Since the kinematics 
around mml appear to be dominated by the jet, Keplerian sig- 
natures cannot be expected for that subregion. The situation is 



less obvious for mm2 (Figure 10 1 where the velocity structure 



at least does not disagree with Keplerian rotation. Regarding the 
alignment of axis, the proposed rotational axis of the large-scale 
toroid (Sandel let al.[ |20"03"] l and the structure around mm2 are 
approximately aligned whereas the jet-axis dominating mml is 
almost perpendicular to that. Unfortunately, such low-number 
statistics do not allow us to derive further conclusions from that. 



5. Conclusions 

Very-high-resolution mm continuum and spectral line observa- 
tions of the two high-mass disk candidates NGC7538IRS1 and 
NGC7538S reveal intriguing information about the small-scale 
morphology and kinematics of these two regions. 

NGC7538S appears as a relatively typical source that frag- 
ments down to the smallest resolvable scales. The large-scale 
single-dish gas clump forms an elongated torus of ~30000 AU 
( |Sandell et al.||2003] ) that fragments into three cores with separa- 
tions on the order of 5000 AU. At even higher spatial resolution, 
these cores show additional substructure and the most massive 
one fragments even further. These data are consistent with hi- 
erarchical fragmentation. While the kinematics of the main mm 
peak mml appears to be strongly influenced by the jet/outflow 
emanating from the source, a spectrum extracted toward the cen- 
tral peak mmla is indicative of additional unresolved rotational 
motions. Higher-resolution data are needed to resolve that. The 
spectral lines toward mm2 also exhibit a velocity gradient, and 
although barely resolved, the data are consistent with Keplerian 
rotation around a low- to intermediate-mass object. Therefore, 
in NGC7538S we are witnessing the formation of a very young 



cluster where the sources within mml have the potential to form 
a high-mass star at the end of the evolution. An additional in- 
teresting feature is that mml and mm2 are strong spectral line 
emitters whereas mm3 is not. While mm2 and mm3 appear very 
similar in the continuum emission, the strong diversity in the 
spectral lines indicate different evolutionary stages. Hence, even 
within areas of ~ 10000 AU diameter, we find cores that are likely 
not evolving coeval. Determining physical reasons for such dis- 
crepancies is beyond the scope of this paper. 

NGC7538IRS1 remains a single source even at -800 AU 
spatial resolution (~ 0.3")- This is even more surprising con- 
sidering that the source is embedded in an already existing near- 
infrared cluster. NGC7538IRS1 has extremely large gas and dust 
column densities corresponding to visual extinction values on 
the order of ~ 10 5 mag. The fact that we still see the central 
source in the infrared implies that the jet/outflow from that re- 
gion should be aligned closely to the line of sight allowing us to 
glimpse through the outflow cavity close onto the central source. 
The central 2000 AU around the source contain a large gas mass 
on the order of 50 M Q , implying central average densities in the 



regime 



of 10 9 . 



Since the position-velocity diagrams of NGC7538IRS1 are 
distorted by the absorption, interpretation of kinematic signa- 
tures are more difficult. Nevertheless, we clearly identify a ve- 
locity gradient in northeast-southwestern direction, consistent 
with the proposed mid-infrared disk emission orientation ( |De| 
Buizer & Minier 2005 ) and perpendicular to the outflow axis. 
Our data do not support the proposed rotational axis based on 
CH 3 OH maser emission ( |Pestalozzi et al.| |2004| |2009| l that is 
inclined to our observed axis by approximately 60 deg. 

At ~0.3" spatial resolution, almost all observed spectral lines 
reveal strong absorption signatures toward the peak of the mm 
continuum emission (that coincides within the errors with IRS1) 
in NGC7538IRS1. While some lines, in particular the lower ex- 
citation temperature lines like those of H2CO appear to be dom- 
inated by blue-shifted absorption indicative of outflowing gas, 
the higher-excitation and higher-density lines exhibit clear red- 
shifted absorption that has to be due to infalling gas. Since the 
jet/outflow is supposed to be aligned along the line-of-sight, it is 
no surprise that infalling and outflowing gas are observed at the 
same spatial position. Estimated mass infall rates are very high, 
on the order of 10~ 2 M Q yr~' . Although we cannot proof that the 
gas will continue to be accreted by the central star, neverthe- 
less, the conditions are sufficient to allow accretion still during 
that already luminous and active phase of the protostellar evolu- 
tion. Combining the large infall rates with the fact of barely any 
fragmentation of the gas and dust core, these data are consistent 
with high-mass star formation proceeding in a scaled-up version 
of low -mass star formation. 

While the presented data already reveal many new insights 
for both regions, significant information is still missing. In par- 
ticular, the proposed accretion disk signatures for both sources - 
NGC7538IRS1 and NGC7538S mml - are "contaminated" by 
absorption and jet signatures, respectively. To overcome these 
issues, one likely needs to resort to even higher excited lines that 
are neither absorbed by the envelope nor emitted by the outflow- 
ing gas. 
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